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Abstract 



Multiple hard and semi-hard parton scatterings in high-energy p + A collisions involve multi-parton correlation 
inside the projectile in both momentum and flavor which will lead to modification of the final hadron spectra relative 
to that in p + p collisions. Such modification of the final hadron transverse momentum spectra in p + A collisions is 
studied within HIJING 2.1 Monte Carlo model which includes nuclear shadowing of the initial parton distributions 
and transverse momentum broadening. Multi-parton flavor and momentum correlation inside the projectile are incor- 
porated through flavor and momentum conservation which are shown to modify the flavor content and momentum 
spectra of final partons and most importantly lead to suppression of large pj hadron spectra in p + A collisions at both 
RHIC and LHC energies. 

Keywords: Cronin effect, HIJING, valence quark number conservation 



In heavy -ion collisions, properties of quark-gluon plasma (QGP) can be studied via jet quenching or suppression 
of high pj hadrons[l] due to parton energy loss[2[[3]. However, initial multiple parton scatterings in cold nuclei can 
also lead to nuclear modification of the final hadron spectra at high-p^. Besides the parton shadowing or nuclear 
modification of parton distributions inside large nuclei, multiple parton correlations inside the projectile nuclei is also 
an important cold nuclear effect. Energy-momentum and valence quark number conservation alone could modify the 
momentum and flavor dependence of final-state parton and hadron spectra at large pj. Hadronization of multiple-jet 
systems from multiple scattering in p + A collisions could also affect final hadron spectra at intermediate and large 
Pt- These cold nuclear effects need to be understood for the study of QGP properties through jet quenching. 

Recently, we find that parton shadowing, initial and final-state transverse momentum broadening, final-state parton 
flavor composition and hadronization of multiple jets, all contribute to the nuclear modification of final hadron spectra 
in p + A collisions|4| within the HIJING Monte-Carlo model||5). HIJING is based on a two-component model for 
hadron production in high-energy p + p, p + A, and A+A collisions. The soft and hard components are separated by a 
cut-off po in the transverse momentum exchange. Hard parton scatterings with pj > po are assumed to be described 
by the perturbative QCD (pQCD), while soft interactions are approximated by string excitations with an effective cross 
section cr soft . Nuclear modification of initial parton distributions should be considered in p + A or A + A collisions. 
HIJING2.0|6 | employes the factorized form of parton distributions in nucleiJT) 



where R^{x, Q 2 , b) is the impact-parameter-dependent nuclear modification factor. 

Besides parton shadowing, the Cronin effect[8| or the enhancement of intermediate pj hadron spectra in p + A 
collisions due to multiple parton scattering should also be considered in our calculation. Multiple scattering inside a 



fa/A(x, Q\ b) = AR A a (x, Q 2 , b)f a/A (x, Q 2 ) 



(1) 
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Figure 1 : (color online) Nuclear modification factor for parton pj spectra in d + A collisions at 200 GeV and p + Pb collisions at 4.4 TeV from 
HIJING2.1 model. 



nucleus can lead to transverse momentum (kj) broadening of both initial and final state partons. This Ay-broadening 
is regarded as the reason for the Cronin effect. In HIJING 2.1, we introduce a Ay-kick in each binary nucleon-nucleon 
scattering to both the initial partons of a hard scattering as well as the final-state partons. The Ay-kick for each 
scattering follows a Gaussian distribution and an energy dependence of the width, 

{k 2 T ) = 0.141og( Vs/GeV) - 0.43GeV 2 . (2) 

After the Ay-kick, longitudinal momenta of partons are reshuffled pair-wise between projectile and target partons to 
ensure four-momentum conservation. The introduced Ay-kick will also influence slightly hadron rapidity distribution. 
After re-tuning the the gluon shadowing parameter s g Q, HIJING 2.1 can describe the rapidity distribution dN c i,/dri 
for d + A collisions at a/* = 200 GeV. We also extrapolate our calculation to LHC energy and give the prediction of 
dNch/drj for p + Pb at -\fs = 4.4 TeV[4|. The experiment data on the charged hadron rapidity distribution at LHC can 
give us further constraints on the gluon shadowing. 

The hard processes in HIJING are described by lowest order pQCD. In p + A collisions, the single jet inclusive 
cross section is proportional to nuclear parton distributions f a /A(x2, p\, b), 



dyid 2 p 7 



K J dy 2 d 2 bt A (b) ^ x x f aip (xi , p 2 T ) x x 2 f a / A (x 2 , pi, b) Ao ~ a ^" :d , (3) 



where, x\j, = pT(e ±y ' + e ±yi )l V^s are the fractional momenta of the initial partons, are the rapidities of the 
final parton jets and the higher order corrections are absorbed into the K factor. The nuclear thickness function is 
normalized to Jd 2 bt A (b) = 1, To study parton shadowing and other cold nuclear effects, we define the nuclear 
modification factor for final-state parton and hadron pj spectra as, 



dN pA ldyd 2 p T 
(N bin )dN pp /dyd 2 p T ' 



RpAtPT) = /A , ['Z ~h±kz > (4) 



where (At) is the average number of binary nucleon-nucleon collisions in p + A collisions. 



In Figjl), we show the modification factor for final parton distributions R p ™ on for minimum-biased d+Au and p + 



J (AW 

Pb collisions at — 200 GeV and 4.4 TeV, respectively, from the HIJING 2.1. We can see clearly the enhancement 
of parton spectra at intermediate pj region. In the default HIJING set, p + A and A+A collision are decomposed into 
independent nucleon-nucleon collisons. Within each nucleon-nucleon collision, hard parton scatterings are simulated 
first, followed by soft parton interactions. In order to maintain the energy-momentum conservation, HIJING subtracts 
the energy-momentum transfers in the previous hard and soft interactions from the projectile. This subtraction will 
restrict the energy available for subsequent binary nucleon-nucleon collisions. Since the time scale of hard scatterings 
is much shorter than that of soft interactions, such a coupling between hard and soft interactions within each binary 
collision might not be physical. In a new set of the HIJING2.1, we turn off this coupling between soft and hard 
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Figure 2: (color online) Quark to gluon ratio as a function of pj in p + p and p(d) + A collisions at 200 GeV and 4.4 TeV from HIJING 2. 1 model. 



parton scatterings. We first simulate all the hard parton interactions in one p + A event, then carry out all the soft 
interactions afterwards. As a consequence, the energy available for each hard scattering is not restricted by previous 
soft interactions anymore. This set is denoted as DHC (de-coherent hard scattering). In Fig|TJ, we can see both the 
nuclear shadowing and the hard-soft coupling can suppress the pj spectra of produced partons. These features of 
parton spectra will be translated into the final hadron spectra after the hadronization process. 

In p + A collisions, the projectile proton will suffer multiple scatterings with the target nucleus. For each binary 
nucleon-nucleon collision there is a finite probability of hard parton scattering, involving independent initial partons 
from projectile and target nucleons. While momentum correlation of multiple partons inside the projectile can be 
neglected beyond the conservation of the total momentum, flavor conservation will limit the availability of valence 
quarks from projectile for each of these hard interactions. This valence quark number conservation will affect the 
flavor composition of produced partons per average binary nucleon-nucleon collision. Shown in Fig|2) are the ratio of 
produced quark pj spectra over gluon, for d + Au collisions at = 200 GeV and p + Pb collisions at = 4.4 TeV, 
respectively, comparing with p + p collisions. We can see that the flavor conservation in p + A collisions can indeed 
suppress the fraction of quark in the produced partons relative to p + p collisions, especially at high pj. Because the 
gluon fragmentation functions are softer than quarks, the increased fraction of produced gluon jets in p + A collisions 
can lead to suppression of final hadron spectra at high pj. 

Besides the momentum and flavor correlation of partons, jet hadronization or fragmentation process can also 
influence the final hadron spectra in p + A collisions. In HIJING default set, the scattered jet partons can form 
jet-shower after initial and final-state radiations, and these radiated soft gluons are ordered in rapidity. While the 
q - q jets shower can form independent string, the gluon jet shower are always connected to the valence quark and 
di-quark of the corresponding projectile and target nucleons as kinks to form string systems. Within Lund string 
fragmentation model[9], all strings are fragmented into final hadron. In p + A collisions, the projectile can undergo 
multiple scatterings, its string system can have more gluon attached comparing to that in p + p collisions. Hadrons 
fragmented from such string system are softer than hadrons from independent fragmentation of individual gluons. To 
illustrate this effect, we also compare the results of independent fragmentation for both p + p and p + A collisions with 
HIJING 2.1. 

In Figj3j, we show the nuclear modification factors for charged hadrons in d + Au at — 200 GeV and p + Pb 
collisions at yjs — 4.4 TeV, respectively, as compared to the existing data from the PHENIX|10| and STAR|11| 
experiment at RHIC. If there is no nuclear effect, R pA » 1 at intermediate and large p T region. The enhancement of 
hadron spectra at intermediate pj is due to the kj broadening through multiple scatterings. As one can see from Fig(3), 
the parton shadowing in nuclei, soft-hard coupling and enhanced gluonic jets due to valence quark conservation all 
lead to suppression of charged hadron spectra. Furthermore, the fragmentation of string systems with multiple gluons 
from multiple hard parton scatterings can further suppress the hadron spectra in the high pj region. 

Since parton shadowing will disappear at large p T \ 12 1 due to QCD evolution and hard scatterings are de-coherent 
from soft interactions, the nuclear modification factor at LHC will likely follow the default set and DHC results at 
low pr and approach DHC without shadowing at large pj with possible further modifications due to hadronization of 
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Figure 3: (color online) Nuclear modification factor for charged hadrons in p + A collisions from HIJING 2. 1 model, comparing with the PHENIX 
and STAR experiment data. See detail in text. 



multiple jets. This possibl scenario is shown as the shaded band in Fig|3j for the nuclear modification R p a in p + Pb 
collisions at the LHC energy. 

In summary, we have studied the nuclear modification of hadron spectra in d + Au and p + Pb collisions at 
the RHIC yfs NN = 200 GeV and LHC energy yfs NN = 4.4 TeV within the HIJING2.1 Monte Carlo model. The 
nuclear modification of the pj spectra is found to be subjected to correlated hadronization of multiple jets and flavor 
conservation in multiple parton scattering in the HIJING model. Experimental test of these effects in the proposed 
p + Pb collisions at the LHC is crucial to disentangle these cold nuclear effects from that caused by jet quenching in 
the hot quark-gluon plasma. 
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